Optically induced coherent intra— band dynamics in disordered semiconductors 
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On the basis of a tight-binding model for a strongly disordered semiconductor with correlated 
conduction- and valence band disorder a new coherent dynamical intra-band effect is analyzed. For 
systems that are excited by two, specially designed ultrashort light-pulse sequences delayed by r 
relatively to each other echo-like phenomena are predicted to occur. In addition to the inter-band 
photon echo which shows up at exactly t = 2t relative to the first pulse, the system responds with 
two spontaneous intra-band current pulses preceding and following the appearance of the photon 
echo. The temporal splitting depends on the electron-hole mass ratio. Calculating the population 
relaxation rate due to Coulomb scattering, it is concluded that the predicted new dynamical effect 
should be experimentally observable in an interacting and strongly disordered system, such as the 
Quantum-Coulomb-Glass. 

PACS numbers: 72.40.+w; 78.47.+p; 72.80.Ng 



I. INTRODUCTION 

Echo phenomena, the prototype being the Hahn spin 
echoi for spin-l/2~systems, rely on the generation of a 
coherent ensemble of excitations with a continuous distri- 
bution of frequencies. After pulsed excitation the macro- 
scopic response of the ensemble decays as a consequence 
of destructive interference effects in the continuum of ex- 
cited frequencies. A second, delayed excitation pulse in- 
duces a rephasing of the individual excited species such 
that at twice the delay time the ensemble shows a spon- 
taneous macroscopic response, the echo. A necessary re- 
quirement for the observability of this coherent dynamics 
are sufficiently weak dephasing interactions on the time 
scale of the pulse delay. The microscopic reason for the 
appearance of an echo is that the second pulse causes 
phase conjugation of the coherent excitation generated 
by the first pulse. In close analogy to the spin echo 
also photon echoes have been observed in ensembles of 
two-level absorbers?:. There are also phonon echoes^ and 
temperature echoesA, even classical mechanic ensembles 
of pendulums can show echo phenomenaS^. 

Photon echoes have also been studied in more compli- 
cated systems, such as ordered^ and disordered^ semicon- 
ductors. The optically excited inter-band transitions in 
a semiconductor cannot be considered as an ensemble of 
independent two-level absorbers due to the strong inter- 
action of the electron-hole pairs. Photon echoes there- 
fore may show a decay as a function of the delay time 
due to the Coulomb interactions, due to disorder, and 
due to combined interaction-disorder effect sSiSiSiiS. In 
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new echo phenomenon has been proposed for 



disordered conductors or Anderson insulators. On the 
basis of a non-interacting one-dimensional tight -binding 
band with diagonal disorder filled with a low density of 
carriers the current response to short externally applied 
voltage pulses was calculated. Assuming excitation with 
two short voltage pulses at i = and t = t, where r was 



chosen larger than the typical elastic scattering time Tei , 
it was predicted that the system spontaneously responds 
with a current pulse exactly at time t — 2t. In contrast 
to spin echoes and photon echoes this current echo is not 
related to phase conjugation. 

The influence of Coulomb interactions on the current 
echo was investigated in Ref>i2,. The numerically exact 
calculation for a small tight -binding system showed that 
the current echo should remain visible in the presence 
of the many-body interaction, however so far no experi- 
mental observation has been reported. 

In a series of papers van Driel, Sipe, and coworkers 
have shown that in semiconductors currents can be in- 
duced optically on an ultra-short time scale using a co- 
herent control schemeW-. In addition, a resulting current 
pulse could possibly be observed using THz-detection 
techniques. 

Stimulated by these results, we studied optically in- 
duced current phenomena in disordered model systems-—. 
By solving the equation of motion for the intra-band 
current in a noninteracting tight -binding model we con- 
cluded that on the basis of the coherent control scheme 
it should be indeed possible to generate a current pulse 
also in a strongly disordered semiconductor. The current 
pulse decays due to elastic scattering. In addition, the 
current traces show signatures of Anderson localization 
(in a one-dimensional noninteracting disordered system 
all states are localized). The application of two delayed 
optical pulses, both generating a current pulse, was found 
to result in a sizable echo in the intra-band current, that 
appears at i = 2t. At the same time the inter-band 
polarization shows the conventional photon echo. 

Continuing our earlier work^^ and making our model 
more realistic we were surprised to discover that the opti- 
cally induced intra-band dynamics in a disordered semi- 
conductor with correlated conduction- and valence-band 
disorder shows features that differ profoundly from both 
the photon and the current echo. By allowing the ef- 
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fective masses of the electrons and holes to differ from 
each other, we find that the related spontaneous sig- 
nals appear at different times. While the inter-band 
photon echo always shows up at t = 2r, the intra- 
band signals split into an advanced and a retarded sig- 
nal current pulse. This splitting depends on the mass 
ratio. For the semiconductor model we find that these 
two current pulses appear only if the disorder is corre- 
lated for electrons and holes, which is often at least ap- 
proximately valid in real low-dimensional semiconductor 
nano-systems where disorder may be due to the fluctua- 
tions of the confinement potential. 



We furthermore find that only the first excitation pulse 
a.t t — has to be designed according to the coherent 
control scheme, which is necessary to generate an intra- 
band current. For the second pulse at t = t a normal 
resonant optical excitation is sufficient. This simple ex- 
citation pulse does not generate a current pulse at t = r, 
however, it is still able to initiate the pair of sponta- 
neous signals. These signals are most prominent for high 
excitation density. Since the many-particle Coulomb 
interaction has been ignored in these calculations, one 
might suspect that Coulomb scattering leads to suffi- 
ciently rapid dephasing such that the new dynamics is 
no longer observable. The consideration of the Coulomb 
interaction for the situation at hand requires to calcu- 
late the response at least up to 5th-order in the external 
light field and to treat the interaction consistently with 
all the relevant correlations (for a 3rd-order treatment 
in a disordered semiconductor model see, e.g.fi2iiS*i&). 
Even for a one-dimensional system this task is beyond 
present computer capacity. In order to get a feeling for 
the relevant time scales we consider a much simpler case. 
For a strongly disordered one-dimensional single-band 
model we calculate the population relaxation time due 
to Coulomb scattering by treating the many-particle in- 
teraction in second Born approximation on the basis of 
Hartree-Fock states. The situation envisaged is the so- 
called Quantum-Coulomb-Glass model, one of the most 
challenging problems in modern many-particle physics. 
The results we obtain can be taken as an indication that 
in fact strong disorder leads to a slowing-down of the dy- 
namical Coulomb effects. For the time being this suggests 
that the current echo should be observable in experiment. 



This paper is organized as follows. In section ITTl the 
model is introduced. The relevant observables, the inter- 
and intra-band currents and the polarization, are defined 
in section Unl The equations of motion are given in sec- 
tion llVl and the coherent control excitation scheme is de- 
scribed in section^ Finally, in section IVTl the numerical 
results are presented and discussed in section IVTll In the 
concluding section IVIIII we consider the possible influ- 
ence of the Coulomb many-particle interaction and give 
evidence for the observability of the predicted signals in 
experiments. 



II. THE MODEL 

We consider a one-dimensional two-band tight- 
binding model with nearest neighbor coupling defined 
to model a direct gap semiconductor. The N sites i have 
nearest neighbor separation \R\ and periodic boundary 
conditions are applied. Every site i at position Ri carries 
two energies and e^. They are distributed randomly 
according to a box-shaped distribution of width . In 
the electron-hole picture all electron energies and hole 
energies as well as the couplings are positive quan- 
tities if the energy is taken to be zero in the gap between 
the valence and the conduction band. The Hamiltonian 
matrix for band A is then 



(1) 



For an ordered situation {ef — Eq, = OeV) we have a 
direct semiconductor with cosine bands and a gap in the 
center (fc = 0) of the Brillouin zone. Close to A: = the 
bands are characterized by effective masses mx related 
to the couplings 

mx = —. (2 

2J^|i?|2 

In the disordered case the relevant disorder parameter is 

In our one-dimensional system all single-particle states 
are localized for non-zero W. The disorder is called un- 
correlated, if the site energies are distributed indepen- 
dently from each other in the two bands. It is called 
correlated if 



J" 



(4) 



where (e^'") are the expectation values of the site ener- 
gies. In this work we exclusively treat correlated disor- 
der. It models, e.g., a disorder potential in a semicon- 
ductor heterostructure with effective dimensionality less 
than three, which is produced by local fluctuations in the 
confining potential. 

The total Hamiltonian then reads 



H 

Ho 



Ho + Hi, 



(5) 
(6) 



Hi = -^E{t)-d. (7) 

a^j, axj are electron creation an annihilation operators, 
respectively. The couplings T in the electron picture 
are unchanged for the conduction band, i.e. T[j = T^^-, 
whereas a change of sign is required for the valence band. 



i.e. = 



_'Y'v 



(8) 
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TABLE I: Parameters of the semiconductor model and the 



hght pulse 




1.316 eV 


r 


34 meV 




3.0 A 


\R\ 


20 A 


Disorder r] 


2 


Number of sites 


71 


Light field described in section FVI 

012 

central wavelength of full-gap-pulse 


7r/2 
1.316 eV 



and 



(9) 



r\xi is the inter-band optical dipole matrix element. 



III. THE OBSERVABLES 



where pij is the inter-band coherence related to the inter- 
band polarization (and the inter-band current) and n^^- 
and n^^ are the intra-band coherences [i ^ j) and densi- 
ties (i = i) related to the intra-band current. 

The intra-band current J™*'"''^ the inter-band current 
jmter ^ and the total polarization P are given by 



2eR 

h 



2efr 



h 

2erv 



J=5^Im[<+J 

i 

- J-5]lm[<M] 

i 



(12) 
(13) 



{r + J^)^Im[p,_i,, +p^+l,i] , 



p 



(14) 



The total polarization is given by 



P = 



V 



and the total current by 



ie 



A ij 

-J- E (r^i -f^^{al^ax.,) 

W ij 
AA' i 



AA' 



(10) 



X ((a+aA'»+i) + (a+aA'.j-i)) 



where R — Ri+i—Ri. These observables have both intra- 
and inter-band contributions due to the first and second 
term in Eqn. Q, respectively. 

Denoting the expectation values of the operators by 







(11a) 






(lib) 




= Sij - {a'l^a^i). 


(He) 



respectively. In Eqn. (|14|) the first and second terms refer 
to the intra- and inter-band polarization, respectively. 



IV. EQUATIONS OF MOTION 

Using the Heisenberg equation of motion and taking 
the expectation values everywhere we obtain 



^ J" (Pi-lj +P-*+lj) + J^J" {Pij-l +Pij + l) 



(15) 



+ j^eE{t) (^R^ ' Rj^ Pij + feu {ntj + n''^, - 5,j) 



d h 
-rKi 
dt 



h 

i 

h 

i 

h 

i 

h 



e^)n^ 



(16) 



-eE{t) {R, - i?j>i^j + fcv {prj - P*ji) 



(17) 



E(t) \(r, - i?,) n% + fc, (p*^- - pj,) 



which differ from the conventional optical Bloch equa- 
tions for a noninteracting tight-binding modeUSti^ii^ by 
the terms containing the position Ri of the sites. 

We numerically solve these equations using the stan- 
dard 4th-order-Runge-Kutta algorithm. 
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V. CURRENT EXCITATION BY COHERENT 
CONTROL 

We apply the coherent control scheme developed by the 
Toronto Groupie in order to optically generate a short 
intra-band current pulse in both bands. The first exci- 
tation at time t = is chosen to be due to a light field 



E{t) = Eie V*ii^ cos 

+ E^e^i'ii} cos [cjt + 01^2] , (18) 



where En are the amplitudes, t^n the temporal widths 
and 01,2 is the relative phase of the two contribu- 
tions having frequency lo (called full-gap pulse) and u /2 
(called half-gap pulse) . We take uj to be larger than the 
band gap, whereas uj/2 is smaller than the band gap. 
In a previous papesi^ we have shown, using the above 
model, that this excitation results in a current depend- 
ing on sin[0i_2] that decays due to scattering at the dis- 
order. This decay is modulated by oscillations which are 
a fingerprint of Anderson localization. It was also shown 
there (see Fig. ^ that a second identical pulse, arriving 
at time t = t, leads to a spontaneous current echo at time 
t = 2t having opposite direction. In that work identical 
electron and hole masses have been taken, i.e. J'^ = . 
The process of generating a current is of at least third 
order in the field amplitudes. The current echo requires 
additional excitation with delayed pulses and is therefore 
of higher order. The analysis presented in the following 
section shows that the appearance of an echo is at least 
of fifth-order in the external fields. Although the current 
pulses have been generated optically, the dynamics seems 
to strongly resemble the current echo as it was originally 
suggestedii for a single band model with current gener- 
ation by voltage pulses. As will be shown in the next 
section, however, the intra-band dynamics initiated in a 
semiconductor is profoundly different from the current 
echo in a single band situation. 



VI. NUMERICAL RESULTS 

A. Dependence on disorder and on electron— hole 
mass ratio 

There are three main findings which point out that 
the dynamics initiated by the coherent control excitation 
characterizes a new coherent phenomenon. 

(i) A spontaneous signal pulse is completely absent if 
instead of correlated disorder according to Eqn. (0J we 
consider a model with uncorrelated disorder. It should 
be mentioned, that in this case for a single disorder re- 
alization current fluctuations are excited even by normal 
band-band excitation. These can be suppressed only by 
extensive configurational averaging over a large number 
of disorder realizations. After sufficient averaging no de- 




FIG. 1: The negative peaks are spontaneous responses fol- 
lowing a succession of two positive current pulses generated 
by coherent control. The effective masses of electrons and 
holes are identical. Note the correct shift of the spontaneous 
signal with 2r. The exciting light pulses have the pulse area 
An = 0.647r (section IVI B|l and a duration of tLn ~ 20fs. 
Each curve has been averaged over 150 disordered realiza- 
tions. From Ref.ii . 



tectable signal can be seen in the resulting current traces. 

(m) The spontaneous response results even for a sim- 
plified excitation sequence. The first pulse at t = is 
taken to generate a current, i.e. it is given by Eqn. H18|l . 
However, for the second pulse at i = r we take a simple 
full-gap pulse with only a single central frequency lo. In 
order to have an influence on the dynamics of intra-band 
quantities, this pulse has to enter the equations of mo- 
tion at least in second order. As a result we still see a 
spontaneous response at the time t = 2r, see Fig. 
We conclude, therefore, that the spontaneous signal does 
not require excitation by two current pulses. Instead only 
one current pulse is necessary and the second light pulse 
can be a single full-gap pulse. The current pulse has 
to precede the full-gap pulse, since the reverse order of 
pulses does not yield a spontaneous response. Because of 
this finding we avoid to call the spontaneous response an 
echo in the following. 

(Hi) For different electron and hole masses Eqn. ||2J), 
J^/J'^ < 1, we find that, surprisingly, there are two sepa- 
rate intra-band responses, one preceding t = 2t and one 
following t = 2t. In particular, the delayed contribution 
is due to the valence band and appears at 



ty 



J" 



(19) 



while the preceding contribution is due to the conduction 
band and appears at 



J" 



(20) 



Note that the inter-band photon echo in all these cases 
always appears at t = 2t, as shown in Fig. 
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FIG. 3: Dependence of the amplitude of the spontaneous sig- 
nal on the pulse areas. Solid line: variation of 1. half-gap- 
pulse, A2 — O.Itt, A4 — 0.2tt; dashed line: variation of 1. 
full-gap-pulse, Ai — O.Itt, A4, — 0.27r; dotted line: variation 
of 2. full-gap-pulse, Ai = O.Itt, A2 = O.Itt. tLn = Wfs in all 
cases. The data has been extracted from averages over 64 dis- 
order realizations. The effective masses of holes and electrons 
are equal. 



VII. DISCUSSION 



FIG. 2: a Spontaneous response in the intra-band current 
after excitation at t = with coherent control pulses and at 
t = 800fs with a single full-gap pulse. For different electron 
and hole masses the spontaneous response splits into a pair 
of signals according to Eqns. 1191 . 12UII . b The inter-band 
current (and also the polarization) for the same excitation 
scenario continue to rephase at t = 2r. Pulse areas are Ai = 
O.Stt, A2 — O.Stt, A4 = TT, the pulse duration tLn ~ lOfs. The 
curves are averages over 64 disordered model realizations. 



B. Dependence on excitation density 



The amplitude of the spontaneous signal depends on 
excitation intensities of the various pulses in different 
ways. Fig. O shows the amplitude as a function of 
the pulse area An of the pulse No. n, defined by 
An = \eE -fcv/hl J^^dt exp[-{t/tLn)'^]- A pulse with 
An = TT corresponds to complete inversion of a two-level 
absorber excited resonantly by pulse No. n. 



It is seen that in particular the half-gap contribution 
has to be strong enough, while the first full-gap pulse 
does not need to have such high intensity. While in the 
limit of low excitation intensity the amplitude of the sig- 
nal depends linearly on the area of the first full-gap pulse, 
the dependence is quadratic for both the half-gap pulse 
and the second full-gap pulse, reflecting the lowest rele- 
vant order of the various pulses. Therefore the sponta- 
neous intra-band response is at least of fifth order in the 
external light field. 



The appearance of the spontaneous intra-band sig- 
nals in systems with correlated disorder can be easily 
understood on the basis of a simplified model. Let us 
assume that we diagonalize the conduction and valence 
band Hamiltonians, resulting in eigenstates \vi>) and lev) 
having energies e^i, and ecv for the valence and conduc- 
tion band, respectively. The Hamiltonian is then given 
in the Appendix Eqn. IjAip . The optical Bloch equa- 
tions describing the dynamics of the inter-band polar- 
ization p^y/ = (a^^ac^i) and the intra-band variables 
= ^i^iy' - {a-tuO-vu') and nl^, = {a+^acu') after ex- 
citation with field E{t) are given in Eqn. ljA4|l in the 
Appendix. Here we are interested in the response of 
to the excitation sequence. We assume 



nf*/,, and n 



the 



vv' 

and 



have been created by the coherent- 



control pulse and are now subject to the excitation with 
the second full-gap pulse. This excitation is resonant 
in the terms proportional to the inter-band dipole ma- 
trix element /i and offresonant in all terms proportional 
to the intra-band dipole matrix element Di,^' . Thus we 
omit all terms proportional to Dy^i in Eqns. (|21|l and 
H22|l . Note that because of the correlated disorder only 
inter-band dipole matrix elements /2 are nonzero between 
pairs of corresponding states, i.e. there is a strict selec- 
tion rule. For tutorial reasons we consider in the follow- 
ing only two states in each band, i.e. \v'2) and |cl), 
|c2). Fig. 01 The equations of motion for 
iox V —\ ,v' — 2 read 



nZ-y and nl^, 



""■21 I r h> 



(21) 
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FIG. 4: The 4-level-model we are using in section IV ill to 
explain how the second excitation pulse at t = t induces 
spontaneous signals in the intra-band current at later times. 
Note that there is a dipole-moment between the two levels 
vl, cl and the levels v2, c2 only. 
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FIG. 5: a The evolution of phases for the case of microscopic 
polarizations which cause the photon echo, b Phase evolution 
of electron and hole intra-band coherences that give raise to 
the spontaneous signals of the intra-band current. Note, that 
the phases of electron and hole coherences are interchanged 
by the second excitation pulse. 



where = (e^i — ev2)/^ and 5c — (cci — ^c2)l^ with 
<5t),c > 0. These equations of motion describe the intra- 
band dynamics after the second full-gap pulse. Just be- 
fore the arrival of this second pulse the tl^^ and 71^2 have 
acquired phases according to 5^ and (5c, respectively, due 
to their free motion in the interval between the first pulse 
and the second pulse, see Fig. [S] Just after the first ex- 
citation pulse their phases were such that (for the total 



ensemble of M levels in each band) a macroscopic intra- 
band current was present. We will denote these phases as 
initial phases. Now our aim is to find the time when the 
phases of tl^ and "rf^ again equal these initial phases. 
We assume for simplicity that the second pulse has the 
form 8{t — r), i.e. it arrives after a delay time r following 
the first pulse (at t = 0) and is extremely short. From 
the equation of motion for p we find the values of this 
variable at time t = t which enter the driving term on 
the RHS of Eqn. (EH): 



dp 



12 



dt 
dp*2i 
dt 



+ i5c2vi P12 

- iSclv2 P21 



--M-^(t)(n^<+n^<) 
^;S.^(t)«+n^<) (22) 



where Sc2vi = (ec2 + eyi)/h and 6civ2 = (cci + £112)/^- 
If these functions, taken at time t = t, are inserted into 
Eqn. (|21|l . we find that the equation for ni^^" has a driv- 
ing term proportional to —n'^^i''') while that for nf^ has 
a driving term proportional to —ni^iir). Note that due 
to correlated disorder the free dynamics in the upper and 
lower pair of states is identical up to a global scaling fac- 
tor given by 5^ /5c- Consequently, n'^^{t) has acquired 
the initial phase of — "■2i"(t) at time <c = (1 + 6v/Sc)t 
and n2{'{t) has acquired the initial phase of — ?t-i2 (''") at 
time ti, = (1 + Sc/Sy)T. This is schematically shown in 
Fig.EjD. Turning now back to the ensemble of more than 
two eigenstates in the bands we see that all these partic- 
ular terms add up at times ty and tc to an intra-band 
dynamics showing the initially generated intra-band cur- 
rent, while all other terms interfere destructively. This 
explains the new intra-band phenomenon. We also learn 
from this consideration, i.e. the sign of the driving terms, 
why the spontaneous signals have the opposite sign to the 
first current pulse. 

In contrast, the inter-band photon echo relies on the 
phase conjugation of the inter-band polarizations p^^i , 
as shown in Fig. This inter-band-phase conjugation 
leads always to a restoration of the initial phases at time 
t — 2t. If, however, one is interested in the dynamics of 
the intra-band quantities, we have to consider that not 
only the inter-band phase factors are conjugated by the 
second pulse, but in addition also the intra-band phases. 



VIII. CONCLUSION 

The experimental verification or observation of the 
phenomena discussed in this paper should be possi- 
ble in cases where dephasing interactions do not de- 
stroy the electron and hole coherences generated by 
the laser pulses. Interactions with the environment via 
phonons and for elevated excitation intensity, in particu- 
lar. Coulomb-scattering provide the major causes of fast 
loss of phase information in these systems. Therefore 
it is necessary to discuss the circumstances under which 
these dephasing processes may still give enough time for 
the coherent spontaneous signals to occur. 
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FIG. 6: Population relaxation rate, F/A, as a function of 
disorder, W/J, for a single-band tight-binding model with 
= 20 sites. The initial population is a Gaussian centered 
around the lower band edge. A is the mean level spacing and 
J is the nearest neighbor coupling. The curves are parameter- 
ized according to U/J, the Coulomb potential at the nearest 
neighbor distance. 

Coulomb-scattering in ordered or weakly disordered 
systems is a rapid dephasing procesaiii^. The question 
here is, how is this fact modified if disorder gets stronger. 
As we have seen strong disorder proves to be more effec- 
tive for the reaUzation of a current echo. It is expected 
that in the locahzed regime the Coulomb-scattering be- 
tween single particle states will be strongly suppressed 
due to the very small overlap of these states. 

In order to have a feeling of what the interplay be- 
tween strong disorder and Coulomb-scattering might 
be we have performed numerical simulations of a one- 
dimensional one-band quantum Coulomb-glass model^^ 
with half-band-filling. After obtaining the mean-field 
approximation by treating the disorder exactly but the 
Coulomb-interaction on a Hartree-Fock level, we study 
the relaxation process of an initially non-equilibrium oc- 
cupation probability distribution over the Hartree-Fock 
basis towards equilibrium caused by Coulomb-scattering. 
This was treated within the second order Born approxi- 
mation. The initial non-equilibrium distribution meant 
to be 'generated' by a laser pulse. 

Increasing the strength of disorder we see a very fast 
decrease in the relaxation rate as can be seen in Fig. 
This effect is attributed to the fact that with increas- 



ing disorder the single particle states become more local- 
ized and therefore reduce the probability of Coulomb- 
scattering. At the end of the relaxation process the 
equilibrium distribution is achieved to be a Fermi-Dirac 
distribution with an effective temperature that increases 
with the excitation energy towards the band center. 

The fact that the efficiency of Coulomb scattering can 
be reduced due to the presence of disorder, may make 
it possible to have the intra-band dephasing times suf- 
ficiently long in order to observe the spontaneous cur- 
rent response predicted here. Ideally suited experiments 
should be performed at low temperatures to reduce scat- 
tering with phonons on samples that are characterized by 
not to weak disorder which needs to be correlated in the 
valence and conduction bands. Semiconductor quantum 
wells with a significant amount of well width fluctuations 
seem to be good candidates for the experimental verifi- 
cation of the phenomena predicted here. 



APPENDIX A 

The Hamiltonian for a system of M states with ener- 
gies ej^ being occupied in the ground state and M upper 
unoccupied states with energies ej; reads 

M 

Ho=Yl XI^^'^^Al^oa,.- (Al) 

\—c,v u—1 

In the electron-hole picture t"^, and are positive quan- 
tities. In the electron picture we have ej^ = —t^ and 
ej; = Cy. The light-matter interaction is given by 

Hl = -E{t) ■ d (A2) 

with the electric field E{t) and the dipole operator 

M M 

d = fl'^(a^^a.t,i, + H.C.) + ^ ^ Di^^'a^^ax^, , 

(A3) 

where D^,ji is the intra-band dipole matrix element be- 
tween pairs of upper and lower states and ft = —ercv 
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dt h 
dt h 



7=1 

M 

7=1 

M 



dt h^"" 



-ii-E{t){p:,, 



(A4) 



7=1 
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